home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / STRCO.z / STRCO
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. SSSSTTTTRRRRCCCCOOOO((((3333FFFF))))                                                            SSSSTTTTRRRRCCCCOOOO((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      STRCO   - STRCO estimates the condition of a real triangular matrix.
  10.  
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  13.       SUBROUTINE STRCO(T,LDT,N,RCOND,Z,JOB)
  14.  
  15. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  16.      On Entry
  17.  
  18.      TTTT REAL(LDT,N)
  19.         T contains the triangular matrix.  The zero
  20.         elements of the matrix are not referenced, and
  21.         the corresponding elements of the array can be
  22.         used to store other information.
  23.  
  24.      LLLLDDDDTTTT INTEGER
  25.         LDT is the leading dimension of the array T.
  26.  
  27.      NNNN INTEGER
  28.         N is the order of the system.
  29.  
  30.      JJJJOOOOBBBB INTEGER
  31.         = 0         T  is lower triangular.
  32.         = nonzero   T  is upper triangular.  On Return
  33.  
  34.      RRRRCCCCOOOONNNNDDDD REAL
  35.         an estimate of the reciprocal condition of  T .
  36.         For the system  T*X = B , relative perturbations
  37.         in  T  and  B  of size  EPSILON  may cause
  38.         relative perturbations in  X  of size  EPSILON/RCOND .
  39.         If  RCOND  is so small that the logical expression
  40.         1.0 + RCOND .EQ. 1.0
  41.         is true, then  T  may be singular to working
  42.         precision.  In particular,  RCOND  is zero  if
  43.         exact singularity is detected or the estimate
  44.         underflows.
  45.  
  46.      ZZZZ REAL(N)
  47.         a work vector whose contents are usually unimportant.
  48.         If  T  is close to a singular matrix, then  Z  is
  49.         an approximate null vector in the sense that
  50.         NORM(A*Z) = RCOND*NORM(A)*NORM(Z) .  LINPACK.  This version dated
  51.      08/14/78 .  Cleve Moler, University of New Mexico, Argonne National Lab.
  52.      Subroutines and Functions BLAS SAXPY,SSCAL,SASUM Fortran ABS,AMAX1,SIGN
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.